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Abstract 

We discuss the inverse problem of determining the possible presence of an 
(n— l)-dimensional crack E in an n-dimensional body f2 with n 3 when 
the so-called Dirichlet-to-Neumann map is given on the boundary of fl. In 
combination with quantitative unique continuation techniques, an optimal 
single- logarithm stability estimate is proven by using the singular solutions 
method. Our arguments also apply when the Neumann-to-Dirichlet map 
or the local versions of the D-N and the N-D map are available. 
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1 Introduction 

Consider an homogeneous electrically conducting body il C R™ which might 
contain an unknown inaccessible crack represented by an (n — 1)— dimensional 
oricntablc surface with boundary S CC !1. Electrostatic equilibrium can be 
modeled by 

in n \ s , 

= , on either side of E , (1-1) 
on dil . 

Here u denotes the electrostatic potential, ip denotes the prescribed potential 
distribution on the exterior boundary dfl. The Robin type boundary condition 
on £ has to be interpreted as follows. 

Having chosen one arbitrary orientation for the normal unit field v on E we 
distinguish by the ± sign the boundary values (or traces) of u and its derivatives 
on the two sides of £ and we denote by z^ + , v~ the normals to E pointing to the 
+, — side of \ E respectively. The impedance coefficients 7 + ,7~ on the two 
side are assumed to be non-negative. 
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We consider the inverse problem of determining £ from boundary current den- 
sity measurements d v u corresponding to one or more choices of the prescribed 
boundary potential ip. 

As is well-known, since Friedman and Vogelius [5T], at least two measurements 
are necessary and in fact, in the two-dimensional setting, it is by now clear how 
two suitable boundary measurements can be chosen in order to have uniqueness 
and stability, [SJ |8J Ql W\ I3H 133 • See also Bryan and Vogelius [H] for a 
thorough review and bibliography data. 

When the space dimension n is three, or higher, uniqueness with finitely many 
measurements is known in few cases, we recall the paper by DiBenedetto and 
the first author [9]. Instead, uniqueness is known when full boundary data are 
available Eller that is when the Dirichlet to Neumann map A : tp H> d v u\du 
is known. 

The aim of this paper is to continue the study of the n-dimensional problem, n ^ 
3, initiated by DiBenedetto and A. and Eller, treating the stability issue when 
a full set of boundary data are available. We shall prove under some a priori 
regularity assumptions on the crack £, that the crack depends continuously upon 
the Dirichlet to Neumann map with a modulus of continuity of logarithmic type, 
with a single log. 

Note that such a modulus of continuity is in fact optimal in view of the several 
examples in [19]. Our approach is based on the use of singular solutions. This 
method can be traced back to Isakov [21] and it appears also in [30] in Eller's 
uniqueness proof. The use of singular solutions for stability estimates is by now 
well-established, gj [TO1 HU [HI H31 HS1 [HI HE However the crack problem 
at hand displays several additional difficulties which have required a completely 
novel approach at various crucial steps. 

Let us describe here the main steps of our proof with their specific difficulties. 
First step 

For two cracks Si, £2 we consider the corresponding Dirichlet to Neumann maps 
Ai, A2. We shall establish an identity which relates Ai — A2 with integrals on 
Si U £2 involving jumps of the corresponding potentials u\ , 1*2 and their normal 
derivatives ( Theorem 13. 21 below ), 

(Ai - A 2 )u 1 u 2 da = / {u 2 [d Vl u\]i - [u\]\d Vl u 2 )do + (1.2) 

JT, 1 \T, 2 

+ / ([u2hd V2 Ui - u 1 [d„ 2 u 2 }2)da + / ([u 2 d Vl Ui]i - [uid V2 u 2 ) 2 )da, 

here u\ and u 2 are solutions to (|1.1[) when £ = £1 = £2 respectively and [-]i, [-] 2 
denote jumps across £1, £ 2 respectively. We refer to (|3.3[) and (|3.4[) for precise 
definitions and to Section [0J for a proof. 

This identity can be viewed as the analogue for the crack problem of the so-called 
"Alessandrini identity" [3J [21] for the Calderon problem. However, in this case, 
its derivation is somewhat intricate due to the fact that the common domain of 
definition of u\ and u 2 is fi\(£iU£2). Such a set, despite the regularity of £1 and 
£2, might be rather wild, thus integration by parts becomes a delicate matter, 
which involves also a preliminary study of the regularity of the potentials u±, u 2 
and the evaluation of the possible singular behavior of their gradients near the 
crack edges <9£i,<9£2, (see Theorem 13.11 below) . 
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Second step 



We apply the above identity to singular solutions U\(-) = Ri(-,y) and u 2 (-) = 
i?2 (-,w) defined on a larger domain and each having a Green's type singularity 
at points y, w placed outside of tt. 

Looking at the right hand side of the identity (jl.2p we introduce the function 



Note that / is harmonic in Q \ (Si U £2) separately in each variable y,w. 
Moreover / is controlled in terms of Ai — A2 when y, w are outside f2. 
Next, by estimates of propagation of smallness for harmonic functions we are 
able to bound f(y, w) when y = w approaches points of £1 AE2 (the symmetric 
difference) . Here the technical obstruction come from the fact that propagation 
of smallness can be performed only on connected sets, whereas f2 \ (£1 U £2) 
may be not. Moreover, not all points of E1AE2 may be reachable from the 
exterior of Q and the estimates of propagation of smallness require that points 
be reachable in a quantitative form which involves the use of chains of balls 
whose numbers is suitably bounded and their radii have to be bounded from 
below. Such requirements induce the introduction of sets Vi of points which can 
be suitably reached from the exterior of (see definition (|4.9I) ) and an ad-hoc 
definition of a variation of the Hausdorff distance for closed sets which we call 
I— distance. 

The crucial point here is that under the a priori regularity assumptions on 
£1, £2 we can show that the Hausdorff distance is dominated by the respective 
I— distance (see Proposition 14.411 . 

Third step 

We show that as y = w tends to a point of E1AS2 then f(y,y) blows up. 
The combination of such a blow up bound and the estimate of f(y, y) in terms 
of Ai — A2 obtained in the previous step lead to the logarithmic estimate of 
dfl-(o"i, £2) in terms of ||Ai — Aa||. The blow up estimate of this step requires a 
careful investigation of the asymptotic behavior of the singular solutions Ri(-,y) 
as their pole y approaches to the crack £$,« = 1,2 (see Proposition 18.31 and 
Proposition 18. 4[) . 

2 The main results 

2.1 Notation and definitions 

In the sequel, we shall make a repeated use of quantitative notions of smoothness 
for the boundary of the domain f2 and for the crack £. Let us introduce the 
following notation and definitions. 





(R 2 (;y)[d Vl R 1 (;w)} 1 - [R 1 (;y)] 1 d Vl R 2 {;w))da + 



([R2(;w)} 2 d V2 R 1 (;y)-R 1 (;y)[d V2 R 2 (;w)} 2 )da + 



([R 2 (;w)d Vl R 1 (-,y)} 1 - [R 1 {-,y)d, 2 R 2 (-,w)} 2 )da. (1.3) 
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In several places it will be useful to single out one coordinate direction, to this 
purpose, we shall use the following notions for points x 6 R™, x' G R n_1 ,x" G 
R n ~ 2 ,n ^ 3,x = (cc',a;„),x' = (x" , x n -i), x" = (x'", x n _ 2 ), with a/ G R™ -1 , .t" G 
R"~ 2 ,a;'" G R"~ 3 and cc„, x n _i, cc„_ 2 G R. Moreover, given a point x G R™, we 
shall denote with B r (x), B' r {x), B" r (x), (x) the ball in R" , R"~ 1 , R™~ 2 , R"~ 3 
respectively centered in x with radius r. 

Definition 2.1. Let SI be a domain in R™. We say that dVl is of class C ' 1 
with constants tq,M if for any P G 90 there exists a rigid transformation of 
R n under which we have P = and 

[~l B r „ — {x G B r „ : x n > <p(x')}, 

where ip is a C ' 1 function on B' TQ satisfying the following condition ip(0) = 
|Vx'y(0)| = and \\ip\\co,i(B^, ) < Mr®, where we denote 

ii H I, I, . \<p(x) - <p(y)\ 

Definition 2.2. Given a, < a 1, we s/iaZZ say i/iai an hypersurface S is 
of class C 1 '" with constants ro, M > if for any P G 5*. there exists a rigid 
transformation of coordinates under which we have P = and 

S H B ro = {(x', x n ) G B ro : x n = ip(x')} (2.1) 

where 

(p:B' ro C W 1 - 1 -> R (2.2) 
is a C 1,a function satisfying 

\tp(0)\ = \Vtp(0)\ = and \\<p\\ c i, a(B > , «S Mr , (2.3) 
where we denote 

\\<f\\ci^(B' ro ) = IIvIIl-(b; ) + r o||V<^|| L=c(£ Q + (2.4) 

+ ^ sup . 

We introduce some notations that we shall use in the sequel. 
For any < r < and any < i\ < r 2 < ro we shall denote 

S r = {x G £ : dist(.T, 92) > r} , (2.5) 

£ r = {x G R" : dist(x, £) > r} , (2.6) 

= {.t G R" : dist(x, 5S,) < r} , (2.7) 

O r = {.TGO c : dist(x,0) r} , (2.8) 

S rur2 = {x G K" : n dist(o;,<9Q) r 2 } , (2.9) 

r r = {x G R" : dist(a;, O) = r}. (2.10) 
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2.2 The D-N map 

Wc begin by defining the Dirichlet to Neumann map. 

For any ip S Hi(dQ,), the unique weak solution to the mixed Dirichlet-Robin 
type problem 

Au = , in n \ E , 

u = if , in the trace sense on 9f2 , (2-11) 

d l ,±u ± — 7 ± u ± = , on E ± , 

is given as the unique minimizer of the quadratic form 

Qe(«) = / |Vu| 2 + [j+u^+j-u- 2 (2.12) 



among all u G \ E), u\ga = (p. 

We denote by < ■, • > the L 2 (dn) pairing between #5(<90) and i?"*(9f2). 

Definition 2.3. TTie Dirichlet to Neumann map associated to (|2.11j) is the 

operator 

A : Hi(dCl) -> #"2(90) (2.13) 

defined by 

<Aip,r]>= / Vd-Vw+ / 7 + v + u + + j~v~u~ (2.14) 

/or every (p,r] € Hi (dfl) where u is the solution to (|2.11[) and w 6 i? 1 (fi \ E) zs 
smc/i i/iai u|ao = 

Note that, as an immediate consequence, we deduce that 

A : Hi(dCl) -> H~i(dQ) (2.15) 

is selfadjoint. 

2.3 Assumptions and a-priori information 
Assumption on the domain 

Given ro , M, D > constants we assume that Q C K n and 

O is of C ' 1 class with constants Tq,M (2-16) 

such that dVt is connected. Furthermore, fi is such that 

diam(ft) < D (2.17) 

Moreover, we assume that the crack E is contained into a closed connected 
hypersurface r C f2 such that 

T is C 1 '" smooth with constants r , M (2.18) 

and it diffeomorphic to a sphere. We also suppose that 

E within T is of class C 1 '" with constants ro, M. (2-19) 
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Namely, for any Q £ <9E, there exists a rigid transformation of coordinates 
under which we have Q = and 

E n B ra = {(x',x n ) £ B ro : x n = ip(x'), x n -i > ip(x")} (2.20) 

where 

V> : S," C M'"~ 2 -> R (2.21) 

satisfying 

^(0) = |V^(0)| = and ||^|| C i, a < M. (2.22) 
Assumptions on the crack impedances 

Given a positive number 7, the crack impedances 7 + and 7" of the unknown 
crack E are such that 

7 ± G C ' 1 ^) (2.23a) 

and 

< 7 ± (x) < 7 for any x £ E. (2.23b) 

We shall refer to the ro, M, D, 7 along with the space dimension n as to the a 
priori data. 

2.4 The main results 

We start by collecting our main stability results for the unknown crack and the 
unknown impedance by means of the global D-N map. 

Theorem 2.4. Let fi,Ei,E2 be the domain and the cracks satisfying the a- 
priori assumptions stated above. If, given e > 0, we have that the D-N maps Ai 
and A2 corresponding to the cracks Ei and E2 respectively, satisfy 

IIA1-A2II 1 1 (2.24) 

then 

d ff (E x ,E 2 ) <C|Iog(e)|-" (2.25) 
where C, n > are constants depending on the a-priori data only. 
Corollary 2.5. Under the same hypothesis of Theorem \2.4\ we have also that 
sup{\^(Q)-^(P)\ s.t. P £ Ei r ° , Q £ E 2 ro nB 2C | log(E) |-,(P)} < C"|log( £ )p' 
where C ,rf > are constants depending on the a priori data only. 



Proof. The Corollary follows by combining the result in Theorem [231 and quan- 
titative stability estimates for the Cauchy problem. For the details of the proof 
we refer to [351 Theorem 2.31. □ 
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2.5 Variants 



In addition, we now state some variants of Theorem 12.41 basically relying on 
other types of data availability. We shall omit proofs since they require only 
minimal adjustments in comparison to the proof of Theorem 12.41 
We start by defining the local version of the Dirichlet to Neumann map. 

Let us fix an open neighborhood A po = B Po (xq) n dCl for a fixed point xq G dCl 

i 

and a given p Q > . We introduce the trace space Hq (A Po ) as the interpolation 
space [Hq(A Po ),L 2 (A Po )]i, we refer to Chap.l] for further details . The 

i i 
functions in i? 2 (A po ) might be also characterized as the elements in H^(dCl) 

which are identically zero outside A po (see for instance [36] ) . this identification 

_ i 

shall be understood throughout. We denote with H 00 2 (A pa ) its dual space, 
which also can be interpreted as a subspace of H~?(dCl). We continue to use 

the notation < •, • > for the duality pairing between 7J 2 (A po ) and H 00 2 (A po ) 
based on the L 2 scalar product. 

Definition 2.6. We shall define as the local Dirichlet to Neumann map asso- 
ciated to (|2.1ip and A po the operator 

A A ™ : l4(A P0 ) Hj (A po ) (2.26) 

defined again by 

<A A ™</?, ?7>= / Vu-Vv+ / j + u + v + +7~u~u~ (2.27) 
Jn Jt, 

i 

for every ip, ?y G iJ 2 (A po ) where u is the solution to (|2. 1 1[) and v G H (CI \ S) 

is such that v\gn = V- 

We now consider the global Neumann to Dirichlet map and we introduce the 
following space of distributions oH~~z (dCl) = {?/ G H~^(dCl) : < ?/, 1 >= 0} . 

Definition 2.7. We refer to the Neumann to Dirichlet map as to the selfadjoint 
operator 

N : oH~^(dCl) -> H^(dCl) (2.28) 

such that 

<i 1 ,Ni 1 >= [ \Vu\ 2 + [ j+u +2 +j- U - 2 (2.29) 

for any r\ G qH ~i(dCl), where u G i/ 1 (f2 \ S) is the weak solution to the mixed 
Neumann-Robin type problem 

Au = , in Cl \ E , 

d u u = T] , on dCl , (2.30) 

d^u^ — 7 ± m ± = , on E ± . 

//, 7 + = 7 _ = on T, ± , we additionally require in (|2.30[) the normalization 
condition J gn u = 0. 
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We are now finally in position to deal with the local Neumann to Dirichlet map. 
Denote A' pQ = dfl \ A po . Let us consider the following space of distributions 

H-i(A po ) = {r, e H-i(8Sl) :<r,,^>=0V^ H^ (A' po )}. 

Definition 2.8. We shall define as the local Neumann to Dirichlet map asso- 
ciated to A po the operator 

N A ?o : H-^(A po ) -> ( I-5(A PO ))" C Hi(dQ) (2.31) 

such that 

<r],N A "o v >= [ |Vu| 2 + / -f+u^+^f-u- 2 (2.32) 

for any i] £ oH~% (A po ), where u £ H l (fl \ £) is the weak solution to (|2.30[) . 
Again, if j + = 7^ = on , we further impose the condition J" af2 u = 0. 

The first variants of our main result concerns the case when the Neumann to 
Dirichlet map is at our disposal instead. 

Theorem 2.9. Let the hypothesis of Theorem \2.J\ be fulfilled. If, given e > 0, 
we have that the N-D maps N± and N2 corresponding to the cracks £1 and £2 
respectively, satisfy 

MiVi — iV 2 II 1 1 s^e (2.33) 

then 

^(Sx.Ea) <C|Iog(e)|-" (2.34) 
where C, rj > are constants depending on the a-priori data only. 

Finally, we treat the cases when the measurements can be performed only on 
an open, non-empty subset S of d£l. Such an instance leads to the introduction 
of the local D-N map and the local N-D map. 

Theorem 2.10. Let fi,£x,£2 be the domain and the cracks satisfying the a- 
priori assumptions stated above. If, given e > 0, we have that the local D-N 
maps A 1 P0 and A 2 P0 associated to A pa and corresponding to the cracks Ei and 
£2 respectively, satisfy 

,| A A P0 _ A A P0 |, i ^ £ , 2 35 , 

£(i? 2 (A P0 ),ff 00 *(A P0 )) 

then 

djr(£i,£ a ) <C|Iog(e)|-" (2.36) 

where C, 7/ > are constants depending on the a-priori data only. 

Theorem 2.11. Let the hypothesis of Theorem \2.J\ be fulfilled. If, given e > 
0, we have that the local N-D maps N 1 P0 and N 2 P0 associated to A po and 
corresponding to the cracks Ei and £2 respectively, satisfy 

"< PO -^ APO "^-(A PO ),( off -i(A PO ))')^ 
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then 



d B (Ei,E a ) <C|log(e) 



(2.38) 



where C, 7/ > are constants depending on the a-priori data only. 

The proofs of the last two theorems can be achieved by combining the results in 
Theorem 12.41 and in Theorem 12.91 respectively with the arguments in |13| where 
the authors provided a quite general method which allow to obtain an Holder 
type dependence of a global D-N map from a local one in a larger domain (see 
also [H] for related results). 

Of course, a more general portion U of dft could be used in the above theorems 
with local data. However the stability constants shall necessarily depends on 
the inradius of such a portion li. For this reason, there is no loss of generality, 
in formulating the above theorems in terms of the spherical neighborhood A po . 

Remark 2.12. For the sake of brevity we only discuss here the stability issue 
for the n- dimensional case with n 3. However our arguments and our results 
could be adapted to the 2- dimensional setting. 

3 The direct problem 

We begin our analysis of the direct problem by providing two results of regularity 
near the crack for the solution to p. II) near the crack, which arc collected in 
the Theorem below and whose proof will be provided in Section [5] 

Theorem 3.1. Let u be a solution to then there exist constants C > 

and a with < a < 1 depending on the a priori data only such 



Moreover, for any p 6 (0,ro) there exists a constant C p > depending on p and 
on the a priori data only such that 



As next step, in the preliminary direct problem treatment we derive an integra- 
tion by parts formula for solutions to the crack problem at hand (|1.1[) . 
Let r,-, i = 1,2 be two closed connected oricntable hypersurfaces of class C 1,a 
as in Section [SJ Just for simplicity of exposition we assume that they are 
diffcomorphic to a sphere. 

By the Jordan separation theorem, Ti disconnects R™ into two connected com- 
ponents il~,C/j, the first one being bounded and the second one unbounded. 
Being dfl connected and f2 bounded, we have H~ CC and dfl C Ui. We 
denote f2^~ = Ui fl Q. Denote i>i the unit normal on I\ pointing to its exterior 
ilf. Further more, the exterior normal to dfl will be denoted by v e (or simply 



IMIc°-(£) ^ C . 



(3.1) 



IM|ci>°(£ p ) ^ C P ■ 



(3.2) 
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If, in addition, we have At; £ L 2 (£l \ Ti) then also the one-sided normal deriva- 
tives d£.v,d~v,d v<s v are defined in the distributional sense. 
We also define 

[d v v}i=d+v-d-v. (3.4) 
Theorem 3.2. (The integration by parts formula) 

Let m € H 1 (r2 \ Ej) be the solution to the problem with £ = £, i = 1,2. 

Then, the following identity holds 

(Ai - A 2 )uiU2cfcr = / (u2[d vi ui]i - [ui]i9 1/1 U2)dc r + (3.5) 
an JEi\e 2 



(M2<9„ 2 hi - ui[d U2 u 2 ] 2 )da + / ([^S^mi]! - [uid V2 u 2 ] 2 )da. 

Ea\2i JSinS 2 

The proof shall be given in Section [5] 

Remark 3.3. ./Vote i/iai the integral on the left hand side of (|3.5p should be 
properly interpreted as < (Aj — A2)ui,ii2 >• .AZso , if r)i = A^u,-, i = 1,2, w;e 
a/so have that the left hand side can be written as < rj\, (N 2 — Ni)r) 2 > . 



4 Singular solutions 

In this section we shall discuss and state the upper bound and the lower bound 
for the function / introduced in (|1.3|) and we shall obtain our main result as a 
combination of the two latter bounds. 
Wc begin by introducing the so called Robin function. 

Fix such that fl CC ft, we shall denote with R the Robin function (or Green's 
function of third kind) associated to the problem (|1.1[) . 

A x R(x, y) = —6{x — y) , in Q, \ £ , 

d„±R+(;y)- 1 ± (-)R ± (;y) = , on E± , ( 4 .1) 

with y € \ S. 

We shall denote with Ri and R 2 the Robin functions solutions to (|4.1[) when E 
is replaced by Ei and E2 respectively. 
Let us now define, for y, w £ \ E 



S^iy,™) 



[ (R 2 (;w)[d Vl R 1 (;y)} 1 - [R 1 (;y)] 1 d Vl R 2 (-,w))da + 

•/El\E 2 

[R 2 (;w)d l/1 R 1 (;y)] 1 da (4.2) 



SinS 2 



Sv 2 (y,w)= / (R 1 (-,y)[d„ 2 R 2 (-,w)] 2 -[R 2 {-,w)} 2 d U2 R 1 (-,y))da + 

•/E 2 \Ei 



[R 1 (;y)d, 2 R 2 {;w)] 2 da (4.3) 

EiHEa 
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note that clearly we have 

f(y,w) = S Sl (y,w)-S'E 3 (y,w) . (4.4) 
By Theorem 13.21 we have that for every y, w e Cl \ fi 

f(y,w)= [ (Ax -A 2 )R 1 (-,y)R 2 (;w)da . (4.5) 
Jan 

4.1 Upper bound on the function f 

Given A, I > we consider the cone 

Ci = {x = (x',x n ) : < x n < Al , \x'\ + ^ < 1} (4.6) 

and for any orthogonal transformation R and any point z, we denote with 

RCi{z) = RCi + z , (4.7) 

the rotated cone whose basis is centered in z. 

Given 7 : [0, 1] — > U O r a simple arc, we define the following set 

7 l = (J Bj(7(*))Ufla(7(l)) ■ (4.8) 
te[o,i] 

Denoting with P{^ 1 ) the vertex of the cone i?Cz(7(l)) and given < r < ro, we 
set 

Vi = {P( 7 ') : 7(0) € IV, 7 ; n (E x U E 2 ) = 0} • (4-9) 

Lemma 4.1. There exist do,Zo > smc/i i/iat i/ e?H (Ei, E2) *S Z ^ ^0 

then 

Ei u e 2 c . (4.10) 

This Lemma will be proved in Section [JJ 

We shall use a variation of the Hausdorff distance which we call Z— distance. 
Definition 4.2. We define the /—distance di between Ei and E2 as follows 



<i;(Ei, E2) = max < sup dist(x, E2) , sup dist(x, Ei) >. (4-11) 

UeEindVi .Tes 2 nav, J 

Here, sup 2 . eSina y i dist{x, E2) is understood to be i/ Ei (~1 = and analo- 
gously for sup xeS2navi dis£(a;, E x ). 

Lemma 4.3. There exist d\, Zi, C > satisfying < e?i ^ do, < Zi ^ Zo swcZi 
Z/iaZ, £/ djjij^ii E2) > cZi i/iera 

d il (E 1 ,E 2 )^C>0. (4.12) 

See Section [7] for a proof of this Lemma. 
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Proposition 4.4. Let f2,Ei,E2 be the domain and the cracks satisfying the 
a-priori assumptions stated above and let l\ > be the quantity introduced in 
Lemma 1 6'. ffl Then, there exists a constant C\ > such that 

dH^EaKCidh^Ea). (4.13) 

This is an immediate consequence of the above two Lemmas, details can be 
found in Section [7] 

With no loss of generality, we can assume that there exists a point O G Ei HdV^ 
such that 

d h =d il (S 1 ,E 2 ) = dist(C»,E 2 ) . (4.14) 

Proposition 4.5. Lei f2 6e i/ie set in 1" satisfying the a-priori assumptions 
stated above. Let l\ > be the parameter introduced before and let Q = P(7~2~) 6 

dVi^ be the vertex of the cone -RC\ (7(1)) for a given simple arc 7. Let y = 

2 2 

Q + hv, where v is the RC ^(7(1)) cone axis unit vector. 

2 

//, given e > 0, we have 

ll Al - A2 Vi W -i (a n))^ e ( 4 - 15 ) 
then for every < h < h, we have that 

C 1 h F 

\f(y,y)\^c^- (4.16) 

where < B < n — 2 and h,C,C , F > are constants depending on the a-priori 
data only. 

Also, the proof of the above Proposition is postponed to Section [JJ 
4.2 Lower bound on the function f 

Let us consider O G Ei r\dVi x the point in (|4.14l) . We introduce a point O' G Ei 
which is defined as follows by distinguishing two cases. 

• If = O G Ei is such that dist((3,<9Ei) < -j-, then we consider a point 
0' G Ei so that dist(0,C) = It follows that dist(C, dE x ) ^ %- and 
dist(0',E 2 ) ^ 

• If O G Ei is such that dist(0, <9Ei) ^ %■ then we set O' = O. 

Proposition 4.6. Let CI be a domain in M. n satisfying the a-priori assump- 
tions. Let Ei,E 2 be two cracks in fl verifying the a-priori assumptions and 
y = hv\ {O'). Then for every h, < h < min{f~o, r, h, ^di ± , ^rd^ } we have that 

\f{y,y)\> c^ ~ c^-h^ (4.17) 

where c = min{i, iD ^ Ira },P = an d c 0j Ci, c 2 > are constants depend- 

ing on the a-priori data only. 

The proof is deferred to Section [7J 
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4.3 Proof of the main Theorem 



Wc now give the proof of our main result. 

Proof of Theorem 12.41 By Proposition 18.31 and Proposition 14.61 we have, up 
to a possible replacing of the constant C in (|4.16[) . that 



,C'h f 



Slh 2 - 71 - — \d p 



c 4 



|3-2n 



c 



where C4 = C2 + C3. By choosing h = qd^™ 2 we have that 



(4.18) 



(4.19) 



with q 



1 



2n-3 

Tl-2 



( 2 C c 4 ) 



and C5 = min{i, 5^}. Hence by 



combining the last two inequalities we obtain 



from which follows that 



C'h F ^ 



> c 5 , 



log(c 5 ) 



log(e) 

Finally, by our choice of h we can conclude that 

d h <C8|tog(e)|-" , 

with c 6 = (g- p C"- 1 |log(c B )|)'?^3 and -q = Fp ^ n 2 _ 3) . 
The thesis follows by Proposition 14.41 with C — cqC\. 



(4.20) 
(4.21) 

(4.22) 
□ 



5 Proof of the regularity estimate 

In this section we shall give the proof of the regularity property of the solution 
u to (|1.1[) and its first order derivatives near the crack. 

Proof of Theorem 13. 1L By the arguments in [Ml Chap. 3], u is Holdcr 
continuous with its first order derivatives up to E except possibly at points of 
<9E. The proof is based on the Moser iteration techniques (see for instance [22l 
Chap. 8]) and by well known regularity bounds for the Neumann problem [5J 
p. 667]. 

We now investigate the behavior of u near the crack edge <9£. Fix xq € <9£, 
up to a translation we may assume that xq = 0. Let us consider the following 
change of variables y = Q(x) 

( y" = z" , 

< y„-i = x n -i - tfj{x") , 

[ Vn = Xn ~ <fi(x", X n -1 - 1p(x")) , 

where ip are the C ' a functions introduced in previous section satisfying (|2.2|) - 
and (i2~20| - (j2~22|l respectively. The map <l 6 C^ a {B^ (0), R n ) and its 

inverse E C ,1 ' a (B ro (0),M n ). 
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With respect to the new variables the crack coincides, within Bjs_(0), with the 

4M 

half hyperplane {y n = 0,y n -i < 0}. 
Denoting with 

A(y) = |det J D$- 1 (j/)|(£>$)($- 1 ( 2 /))(£)$) T ($- 1 ( 2 /)) , (5.1) 

i + {y) = 7 + (^ 1 (y)) , T(y) = 7"(^ 1 (y)) , (5.2) 

v(y) = u($- 1 (y)) (5.3) 
we have that v € i/ 1 (BT a _(0)) is a weak solution to the problem 



f div(A(y)V«(y)) = , inB ^(°) 
i i(y)V W ± (y)^ ± - 7 ± (y) U ± (y)=0, on 57^(0) 

^ 4M 

where 

%(0) = J Bn,(0)\{ y „ = 0} 

and 

B7^(0) = S Z!L (0)n{y„_ 1 <0} 

and v+ = (0, . . . , 0, -1), I/" = (0, . . . , 0, 1). 

We introduce the following system of variables z = $f(y) 



(5.4) 



z" = v" , r = Jy 2 n _ 1 + y 2 n 



. /r(r + 2/ n _i) 



z„ = sign{y n )y - , 

For the reader's convenience we express both systems of variables in cylindrical 
coordinates also 

y = (y",r cos 9,r sin 9) (5.5) 

and 

9 9 

z = (y",rcos-,rsin-) (5.6) 

withO <r<^,-7r<0<7r,-a#<^<^#,i = l...,n-2. 

The underlying idea here relies on mapping through ^ the set B ro v^(0) into 

SM 

the half ball i?+^(0) = {z e B ro ^ (0) : z„_i > 0} so that the two side of 
the -B~^:(0) are mapped into the flat part of B~^ o ^(0) (see also [9l Remark 

SM IF" 

C.3.1.]). 

Moreover it can be verified that the map * € V^ 1 '°°(S Ia v™(0),E n ) and vf -1 € 

8 A-/ 

W 1 ' 00 ^ v5r(°)' R ")' 

Setting 



B{y) = \&etDi>- 1 (z)\(Di>){i>- 1 {z))A(i>~ 1 (z))(Di>) T (i>- 1 {z)) , 
io(z) = u(* _1 (z)) , 
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we have that w £ H 1 (B J ^ ^-(O)) is a weak solution to the problem 

SM 

div{B(z)\7w(z)) = , in B + ~ (0) , 

"fer (5.7) 
B(z)X7w(z) ■ v - j(z)w(z) = 0, on B TJi ^_ (0) n {z„_i = 0} , 

SM 

where z>= (0, ...,0,-1,0). 

Observing that P G ^(5+^(0)), 7 G ^(B^^nf^-i = 0}) and dealing 



again as in [Ml Chap. 3] we infer that w £ C°> a (B+ oVS (0)). 

8M 

Finally, coming back to the former system of coordinates, we obtain the thesis. 
□ 



6 Integration by parts, proofs 

In this Section we shall deal with the proof of our "Alessandrini identity" type 
formula tuned for the crack problem at hand. 

Lemma 6.1. (The divergence formula over 17 \ (I"i U I^)) 

Let F be a vector field such that F G C 1 (tt\ (Ti Ul^)) and moreover F G C(fij_) 

and F G C(fi!_) with i = 1,2, £/ien £/ie following holds 

divFdx = / [F ■ u^daix) + [F ■ v 2 ] 2 da{x) + 

Jn Jt^Tz Jr 2 \ri 

+ / [F ■ y x ) 1 da{x) + I F-v e da{x). (6.1) 
Jrinr 2 Jan 

Proof. Given < p < ro, we have that by the compactness of f2 we can find 
a finite number of points Pj, i = 1 . . . , N such that LifL 1 B p (Pi) covers f2. Let 
us observe that due to the regularity hypothesis made on i — 1,2, we can 
choose p small enough so that on each ball B p (Pi) with i = 1, . , . , N, Ti and I^ 
are separately graphs each with respect to a suitable reference system. 
Let {ai}fL 1 be a smooth partition of unity subordinate to the open covering 
Vjf =1 Bp(Pi), namely we are assuming that 

i) < on < 1 , o>i G C^(B p (Pi)) ,i = 1, . . . ,N ; 

") EjIi"i = lonUil 1 P p (P i ). 

Then, we have 

f divFdx = V / div(aiF)da: . (6.2) 
Jn i=1 JB p {Pi) 

The only interesting cases to consider are when B p (Pi) contains TinT 2 , whereas 
in the other ones the divergence theorem may be applied in a straightforward 
fashion. Let us fix a small aperture 8q and let us distinguish two cases. 
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1. Vxq G Bp(Pi) PlTi the tangent planes of Ti, T2 at xq have an aperture 



2. 3xq G Bp(Pi) fl Ti n T2 and the tangent planes of ri,r2 in Xq form an 
aperture 8 < 8q. 

Case 1. In such a case B p (Pi)\ (Ti Ul^) is composed by finitely many Lipschitz 
domains and the divergence theorem can be used in each component separately. 
Note that the same occurs when Fi n T2 n B p (Pi) = 0. 

Case 2. In this situation, if one chooses #0 and p sufficiently small in terms of 
ro, M one obtain that there exists x G Y\ n T2 PI B p (Pi) such that ri,r2 are 
tangential at x. In this case having chosen p sufficiently small, Ti and I^ are 
simultaneously graphs with respect to the same reference system. Moreover, we 
consider the following three domains. 

a) U = {(x',x n ) G Bp(Pi) : x n > max{(pi(x'), ip 2 (x')}} ; 

b) I = {(x',x n ) G Bp(Pi) : u\u\{tpi(x'),Lp 2 {x')} < x n < max{^i(a/),^ 2 (^)}} ; 

c) L= {(x',x n ) G Bp(Pi) : x n < min{ipi{x'),Lp 2 {x')}} . 

Both U and L are Lipschitz domains. The set / may not be Lipschitz and 
disconnected, but it is a normal domain between Lipschitz graphs. Hence in all 
such sets the divergence theorem holds true. □ 



Proof of Theorem 13.21 There exists a sequence of C^°(R n ) functions 
ip m , m G N, satisfying the following properties. First, $J ip m $J 1, ip m is 

identically equal to 1 in U\_ and ip m is identically equal to zero outside IA™ . 

2 m 

Second, we have that 

\V<p m \^Cm \U\ K (6.3) 

and hence 

\Vip m {x)\ 2 , (6.4) 

where C > is a constant depending on the a-priori data only. 
We notice that 

div(w 2 V3 m Vui) = Vu 2 • Vuiip m + Vu 2 • V(p m Ui. (6.5) 

We observe that Vtti • Vu 2 tp m a.e. in f2 as m 00 and also in L 1 (f2) by 
dominated convergence. On the other hand we have 

/ |Vu 2 ■ V^ttil «S I / \Vip m (x)\A [j ul\Vu 2 {x)\A . (6.6) 
Jn / \ u ±. ) 

By the bound in (|6.4[) and observing that u\\S/u 2 (x)\ 2 G L 1 (f2) we can conclude 
by the absolute continuity of the integral that the right hand side of (|6.6p tends 
to zero asm4 00. 
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Hence we found that 

div(u2(x)Viti(x))rfx = lim / div(u 2 (x)(l — ip m )\7ui(x))dx (6.7) 
n m ^°° Jn 

Using the divergence formula (|6.1|) with F = 7/2(1 — ip m )\7ui we get that 



/ div(u 2 (l - Lp m )Vui)dx = I [u 2 {l - tp m )d vi ui]ida(x) + 
Jn Jri\r 2 

+ / [u 2 (l - tp m )d V2 ui} 2 d(j(x) + I [u 2 (l - ip m )d Ul u 1 ] 1 da(x) + 
Jr 2 \ri Jr 1 nr 2 

"2(1 - ¥ra)d Vc Uid(j{x). (6.8) 

an 

Let 5* be any of the portions Ti \ T2, T% \ T\, T\ n T 2 . We claim that 

9m[U'2d Vi ui]idcr(x) -> , (6.9) 

s 

with i = 2 if S = T 2 \ T\ and with i = 1 in the remaining cases. 
We observe that in order to prove our claim it sufficient to establish that | Vui | £ 
i 1 (S'). If <9£i n S — the integrability of |Vui| over S easily follows from 
Theorem 13. II Let us then analyze the case when <9Ei n S ^ and distinguish 
two situations. 

i) 9Ei intersects S trans versally; 

ii) <9£i intersects S tangcntially. 

We begin by analyzing the case i) and observing that in such a case the intersec- 
tion V = 9Ei n S is a (n — 3)— manifold. We can find a finite number of points 
Pi, i = 1, . . . , N in V such that U^L 1 B?(Pi) covers V, where f will be fixed later 
on. After a translation we may assume that P, = and fixing local coordinates, 
we can represent S as a graph of a C 1,a function ip satisfying (|2.1j) - (|2.4j) . Let 
$ G C 1 ^ a (Bss T (0),R n ) be the map defined as follows 

$(y',yn) = (y',y n + <p(y')) ■ (6.10) 

we have that there exists Q\ , 9 2 , B\ > 1 > 9 2 > constants depending on ro and 
M only such that for any r £ (0, jjj) it follows 

S e2r (0) c $(B r (0)) C B SlP (0) . (6.11) 

The inverse map -I^ 1 e C 1,a {B ro (0), M n ) and it is defined by 

x n ) = (x\ x n - (fi(x')) . (6.12) 

Moreover, by our assumptions on <9Ei and by the implicit function theorem we 
have that there exists f > depending on the a-priori data only such that for 
any r £ (0, f) 

$" 1 ( J B e2r (0) HV) C {y £ B' r (0) : y„_i = ^i(yi, . . . ,y„- 2 ), J/n-2 = Mvu---,Vn 
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where V'i £ C^ a (B'J!(0), R) and V2 e C 1,a (i?"'(0), R). In particular, when n = 3 

the set $ _1 (i?e 2r ,(o) D V) reduces to a single point. 

Let * € C 1,a (B'r ,R n_1 ) be the map defined as follows 

*(z"', Z n - 2 ,Z„-l) = {z"\ Z n - 2 + M*'"): Zn-X + M*")) ■ ( 6 - 13 ) 

As before it can be proved that there exist constants 83 , 64 such that #3 > 1 > 
64 > depending on ro and M only such that for any p S (0, jtj) it follows 
that 

B' eip (0) c *(B;(0)) C B' e3p (0) . (6.14) 
The inverse map ty' 1 e C 1 '"^, R" -1 ) and it is defined by 

*-V,2/„- 2 ,2/„-i) = (y"',y n - 2 -My'"),yn-i-My")) • (6.15) 

Let x be a point in Be^^O), then by Theorem 13.11 we may infer that there 
exists a constant Ci > depending on the a-priori data only such that for any 
x eV (~) Be^r (0) we have 

8M 



x — a; 



C|x - i^ 1 sC Cidist(a;, V n B^ (O)) ■ (6.16) 



Let then a; be a point in £? 84? (0) n S and let z' £ B'^_ (0) be such that x = 

8M 4M 

^(^(z'), 0). Furthermore, let x$ € V he such that \x — ajo | = dist(x, V) with 
xo = ^(^(z'", 0, 0), 0). We have that there exists constant C2 > depending 
on the a-priori data only such that 

|V«i($(*(z'), o))| < c 2 |$(*( 2 '), o) - $(*(*"', 0, 0), o)!"" 1 . 

Finally, by the C 1 '" regularity of $ _1 and ^P -1 we can infer that there exists a 
constant C3 > depending on the a-priori data only such that 

|VM<W),0))| «S C 3 |(z',0) - (z'^O)!"" 1 . (6.17) 

From the above estimate we deduce that there exists a constant C4 > depend- 
ing on the a-priori data only such that 

/ Vui($(*(z'),0))|dz' < C4 ■ (6.18) 

JB'^_ (0) 

Hence choosing f = ^2^47^7 and by a covering argument, we obtain that 

/ \Vui(x)\dx < 00 . (6.19) 
Js 

We now treat the case ii). Since in this case the intersection <9£i P\S might be an 
irregular set, we find convenient to consider the orthogonal projection operator 
II : dTi — > S and we define W = II(<9ri) which is an (n — 2)— manifold. As 
before, we can find a finite number of points Pi, i = 1, . . . , N in W such that 
U^L 1 Bf(P — 1) covers W, where f will be chosen later on. Dealing as before we 
can locally flatten the hypersurfacc S by the diffcomorphism <f>. 
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Furthermore, by our hypothesis on 9Si we have that there exists r > depend- 
ing on the a-priori data only such that for any r 6 (0, r) we have 

<S>-\B e AO)nW) c {y' €B' r (0):y n -i = ^(yi,...,y n -2)} , (6.20) 
where if € C^ a (B' 7 ! (0),R). 

Let * G G 1 '° l (B' r ,R n " 1 ) be the map defined as follows 

*(z", *„_!) = (z", z n _! + #(z")) . (6.21) 

There exist constants 8 5 , 9 6 such that 6* 5 > 1 > # 6 > depending on r and M 
only such that for any p G (0, it follows that 

5^(0) c $(S;(0)) c 5^(0) . (6.22) 

The inverse map G C 1 ' Q (B p , R"" 1 ) and it is defined by 

*"V, = (j/",j/„-i -^(y")) • (6-23) 
Let then a; be a point in -Be^O) (~1 S* and let j/ G Be^r (0) n OTi. Arguing as 

8M 8M 

for case i) we can deduce by Theorem 13.11 that there exists a constant C5 > 
depending on the a-priori data only such that 

|Vui(aO|< Calx-yl"- 1 ^ C 5 \x ~ Uiy)^ 1 «C (6.24) 
< C 5 dist(x,W DBe AL (0)) a - 1 

8M 

Hence let z' e B^(0) be such that x = $($(z'),0) and let y G W be such 

4M 

that dist(x, W r nBe4£(0)) = |x-y | with y = $(*(z", 0), 0), then we have that 

8M 

there exists a positive constant C@ > depending on the a-priori data only such 
that 

|V Wl ($(*(z'),0))| ^ C 6 \(z',0) - (z^O)!"" 1 . (6.25) 

Hence dealing as for the case i) and fixing the radius f = 82 0$ j^j- we get that 
also in this situation \Vui(x)\ G ^(S). 

Hence, our claim (|6.9|) is proved. Combining (16.7[) . (|6.8[) and (|6.9[) we get that 



/ (9 I / c uim 2 - d„ s u 2 ui)da = / (u 2 [9 l / 1 ui]i - [ui]i9 I y 1 U2)rfcr + (6.26) 
Jan Jri\r 2 



+ / (M2<% 2 Mi - u 1 [d^ 2 u 2 }2)da + / ([^S^MiJi - [uid V2 U2\2)d(J. 
JT 2 \Ti JT 1 r\T 2 

Noticing that the integrals over 1^ \ Sj, i = 1, 2 cancel each other since [ujjj 
and [<9u;]i vanish there, the formula (|6.26|) can be simplified as follows 

/ {d Vts u\U2 — d VB U2Ui)da = / {u 2 [d Vl ui\i - [ui]id„ 1 U2)da + (6.27) 
Jan JSi\s 2 

+ / ([u2hd V2 ui - u 1 [d„ 2 U2}2)da + / ([u 2 d Vl Ui]i - [uid V2 u 2 ]2)da. 
Je 2 \Si JSins 2 

Finally, the desired identity follows by selfadjointness of the Dirichlet to Neu- 
mann map A : Hi(dfl) -> H~i(dfl). □ 
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7 Proof of Proposition 14.51 



In this section we shall provide the proof of the Proposition 14.51 together with 
the related auxiliary results stated in Section @] 

7.1 The /-distance 

Proof of Lemma 14.11 We give a sketch of the proof based on three steps. 

1. Being Ej, i = 1, 2 contained into a C 1,a hypersurface Ti, i = 1, 2 and by 
the arguments carried over in [Jj Proposition 3.6] we may infer that there 
exist number g?o,Po:^o > 0>0 < Po < r o for which the ratio ^ only 
depend on a and M, such that if we have 

dH(Ei,S 2 )<do (7-1) 
then for any P S Ei we have that 

T t n B po {P) = {xe B pa {P) s.t. x n = ifii{x')} , i = 1,2 (7.2) 

and ||(^i - (p 2 \\ c i,f, B < Cr 2+2Q dp +2Q , where C > depends on a 

and M only. 

Moreover, the functions max{y>i, c/?2} and mhi{ipi, (^2} are Lipschitz with 
Lipschitz constants bounded by L = 2M + CrJ +2a dQ +2a . 

2. We recall that in our regularity hypothesis for any P £ E;, i = 1,2 we can 
define two unit normals Vi{P) and —Vi(P) with i = 1,2 according with 
the criterion stated in Section [21 Moreover, for any point P € E^, i = 1,2 

III 10. 

we can find a set 7j 2 with A = ^7 so that Ej n 7/ =0 and the axis of 

To. 

the cone RCl^ (within 7, 2 ) is oriented along Vi{P), i = 1,2. The same 

2 

holds true for the unit normal — i/i{P) i = 1,2. 

3. Let us assume that dni^i, E2) ^ do and let P 6 Ei U E2, say for instance 
P G Ei. Without loss of generality, we may assume that in the local 
representation of Ei and E2 as relative graphs P belongs to the graph of 
ma,x{ipi, ^2}. We notice that up to replacing ^ with l = min{^-,^-} 
and A with ^ we have that the set 7^° introduced in the previous step is 
such that 7^ n E2 = and this concludes the proof. 

□ 

Proof of Lemma Hm Let di =min{^,^}. Ifdjj(Ei,E 2 ) ^ d x then without 
loss of generality we may infer that there exists Q e E 2 fl [E\ d ). 
We can find a continuous path 7 : [0, 1] — > E\ dl such that 7(0) <E T r and 
7(1) = Q. Let now Q' = 7(f) with t = mf{t : j(t) £ E 2; J and let P' E E 2 such 
that \P' - Q'\ < li with l\ = min{^, 4^}. Furthermore, let us denote with 7' 
the restriction to [0,t\ of 7. Up to a possible replacement of the constant l\ we 
have that P' £ dVi 1 (where the construction of such a Vi 1 is based on the path 
7' and A = 217). Finally we get 

dist(P',Ei) > di -h > (7.3) 
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which implies that 

d h (H 1 ,H a )^di-li>0. (7.4) 
The thesis follows with C = d\ — l\. □ 

Proof of Proposition 14.41 We distinguish two cases 

1. dff(Ei,E 2 ) ^ d , 

2. dij(Ei, E 2 ) > do, 



where do is the constant introduced in Lemma [ 
Case 1. In such a case we have by Lemma [6.91 that Ej n dV^ = Ej, i = 1,2 
and hence dH-(Ei,E 2 ) = <i; 1 (Ei,E2). 

Case 2. In this situation we have that being dff(Ei,E 2 ) > d\, by Lemma l4~3l 
we can infer that d/ 1 (Ei, E 2 ) > C. Hence we get 

daCEx.Sa) < dg(S ^ S2) d il (E 1 ,E 2 ) < ^d il (E 1 ,E 2 ). (7.5) 

where I? is the a- priori bound on the diameter of VL introduced in (|2.17|) . The 
thesis follows by choosing C\ = maxjl, □ 



7.2 Proof of the upper bound on / 

Proof of Proposition 14.51 Let 7 be the simple arc in the definition of Vi± and 

2 

let x\ = 7(0) 6 T r with < r < Let us fix y £ cS|,4 r , where Sr,4 r is the 
set introduced before. Let us consider f(y, ■) and let 7 be the simple arc in the 

definition of Vi^ , then we have that 

2 

A w f(y,w)=0 in Of, (7.6) 

For w £ Sr, 4r , by (j2~2^)) . ([43]) and ([571^ 

|/(fcu;)|<C||Ai-A 3 ||=Ce. (7.7) 

By the three spheres inequality for suprcmum norms of harmonic function we 
have that there exists a constant < r < 1 such that 

ll/foOIU-CB^O) < ll/^OIIL-^^jjll/^OIIi-V^)) ■ (7.8) 

We consider a point zz> lying on the arc 7 and such that zz> £ 7 "2" \ -B 2 /i(<3). 
Let us define {2^}, i = 1, . . . , s as follows, x\ it has been already introduced, 
.?*;+! = 7(ii) where ij = max{t.|7(t) — Xj] = r} if \xi — w\ > r, otherwise let 
i = s and stop the process. By construction, the balls Bz(xi) are pairwisc 
disjoint, \xi+i — Xi\ = r for i = 1, . . . , s — 1, \x s — w\ ^ r. By (|2.17[) we have 
that there exists a positive constant /3 depending on the a priori data only such 
that s < (3. An iterated use of the three spheres inequality for f(y, •) gives that 
for any < p < r we have 

\\f(y,-)h~(B P ( W}} < limOIII-^.^ojH/Cv,-)!! 1 -^^ • (7-9) 
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l 1 

We observe that for w £ 7^" \ i?2ft (Q) we have that 

\S El (y,w)\ 

^ / (\R 2 (x,w)[d Vl Ri(x,y)]i\ + \[R 1 {x,y)] 1 d„ 1 R 2 (x, w)\)da + 
|[i? 2 (x, tw)9 I / 1 i2i(x,y)]i|do-(x) ^ 

Sins 2 

< C[( (\x~w\ 2 - n + \x-w\ 1 - n da(x)+ f \x - w\ 2 - n )da{x)\ < 

< C/i" 1 . (7.10) 
Similarly, we get \Ss 2 (y,w)\ ^ Ch^ 1 . Then we can conclude that 

\f(y,w)\ ^CTT 1 for any w £ \ B 2h (Q) . (7.11) 
Hence, we have that by (|7.7[) and (|7.11l) 

ll/(v,0IU-<B fi («0) ^Ce'V" 1 . (7.12) 
We introduce the following set of quantities for k ^ 2 

= arctan (^~^J > ® l = arct an > tf.ldi) 
1 + sm(fli) 

Ai = min I —7-7^- , ■= \ , A fe = x^k-i (7.15) 
|_ 1 + sm{&) 3sin(0) J 

pi = Ai sin(0i) , p fc = xPfe-i (7-16) 

wi = Q + Ai? , w k = Q + \ k v. (7-17) 

By repeating the argument outlined in [101 Proposition 3.5] (see also [7]) and 

based on iterative application of the three spheres inequality over a chain of 
balls B pi (u>i), . . . , B Pk(r) (u>£,( r )) within the cone we obtain that 

||/(y,-)IU-(B Pfe(r) (^ M )) < Cs^^- 1 ^ 3 - 1 , (7.18) 

where fc(r) is an integer such that k(r) ~ 1 1 j^^j r ^ with < r < minl^j 1 , ro}. 
Let us now consider f(y,w) as a function of First, we observe that 

A y f{y,w) = in Of;, for any w e n|. (7.19) 
For y,w € 7"^ \ Bh{Q), y ^ w we have that 

|5 El (y,«;)| s? / |i?a(a:,«;)[^iJi(a; jy )]i| + 1^1(0;, V )]i^iZ2(a:,ti;)|de7 + 

JEi\S 2 

|[i? 2 (a;,ty)^ 1 ^i(x,2/)]i|d(T(a;) < 

Eins 2 

s$ C [ \x- w\ 2 ~ n \x - y\ 2 ~ n + \x - y\ 2 ~ n \x - w\ 1 ~ n da{x) + 
+C [ \x - w\ 2 ~ n \x - y\ 2 - n da(x) 
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Moreover, dealing as in Proposition 18.31 we get 

in / \i ^ f h~ x \ log \y — nil I if n = 3 , 

r>|, !p-« ifn>3. ^ 2 °) 

and similarly for Ss 2 (y,w). Therefore, 

I/(V,«0I<{ f^'V™' 1 if f " = 3 ,' (7-21) 

with y,w e 7"2" \ Bf l (Q). Moreover, for y e 5r 4 r and w € 7^" \ Bh{Q) using 
([T~fg]) we have 

l/fe.^KC^-V- 1 . (7.22) 

Proceeding as before, let us fix io G 7"^ such that dist(w, Q) — h and y G £5.4^ 
Again, taking y± = Q + \{u and using iteratively the three spheres inequality 
we have 

\\f(;w)h^ B (V1)) < II/C^IILL^^II/C-,^)!! 1 -" 8 . (7.23) 

where r and s arc the numbers established previously. We now distinguish two 

cases 

i) n = 3 

ii) n > 3 

We begin by analyzing the case i). 

By combining (|T2Tj) . (f?T2"2l and (j7T2"3]) we have 

H/O.wJIU-CB^CvO) < Ch r "- 1 \]ogh\ 1 - r ' E - r *T hW - 1 . (7.24) 

We observe that for h sufficiently small we have that | log h\ 1 ~ T ^ h~^ T . And 
hence from the above estimate we deduce that 

\\f(;w)\\ L ~ (BrM) < CT^-V 28 ^ 1 . (7.25) 

2 

Once more, we apply iteratively the three spheres inequality over a chain of 

balls contained in the cone ^^(7(1)) and we obtain 

2 

H/(^)ll*~(^ w <*<*») < Cfc(K--Da-r^-) e r-^-«. (7 . 26) 

From the above inequality, choosing y = w = Q + 2hv we have that 

|/(„)KcrV VWkH1 (7.27) 

where -B = 1 — ^t 2 ". We observe that, for < h < ctq with < c < 1 depending 
on the a-priori data only, we have k(h) «C c\ \ogh\ = —clog h, so we deduce that 

r 2k(h) > exp(-2clog/ilogr) = h F 
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with F = 2c|logr|. 
Finally we obtain that 

\f(y,y)\ < h -B £ r>^-» ^ Ch - B e^- 2hF . (7.28) 

Hence the thesis follows with h = cro, C = r 2/3 ~ 2 . 

For the case ii) the estimate (|7.28[) holds true with B = n — 2 — r^(n — 3) — r 2/3 
the other constants remaining the same and can be achieved by adapting the 
argument above. □ 



8 Proof of Proposition 14.6 




We premise the proof of Proposition 14.61 with several preliminary results. 

Lemma 8.1. There exists a constant C > depending on the a-priori data 
only, such that the weak solution v £ \ E) to the problem 

tnfl\S , 
on dft , (8.1) 
= , on E ± . 

is such that v{x) ^ C in fl. 

Proof. The existence and the uniqueness of the weak solution v £ 1 (f7 \ E) 
to the problem (|8.ip is a consequence of standard theory on the boundary value 
problem for the Laplace equation and the non negativity of the coefficients 7 + 
and 7~ . We understand that v satisfies 

/ Vv ■ Vip dx + / r y + v + p + da + / r y~v~<p~da = / tp da (8-2) 
Jn\s Jt, Jt, Jan 

for any tp £ H 1 ^ \ E). 

Let V- £ i7 1 (£1 \ E) be the negative part of v, namely t>_ = — min{t;,0}. 
Choosing tp = v_ in (|8.2[) we have that 

/ \\7v^ 2 dx+[ j + \v±\ 2 d<T+ j~\vZ\ 2 do- ^ . (8.3) 

Then by the Poincare inequality we deduce that V— = a.e. in £1 and hence 
v(x) ^ a.e. in fi. 

Let G 9(fi \ E) be such that imnv(x) = v(xq). 

As a consequence of the Giraud's maximum principle (see [221 Theorem 5]) we 
have that xq £ E. Without loss of generality, we may assume that if xq £ E 
then mmv(x) = v + (xq). 

Let us denote for any < p < ro with 



Ap(x ) 



B p (x ) n 0+ if x G E 
B p (x ) \E if x G <9E. 
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By the weak Harnack inequality at the boundary (see [35J Lemma 3.2]) and by 
the non negativity of v we have that there exist a radius r, < r < r$ and a 
constant C > depending on the a-priori data only such that for any < p < r 
we have 

v(x )= min v{x)^ C\\v\\ L 2 (A {xo)) . (8.4) 

x£A p (x ) 

Moreover, dealing as in the proof of Lemma 3.3 of [35] and relying on an iter- 
ated use of the Harnack inequality wc can conclude that there exists a constant 
C > depending on the a priori data only such that v{xq) ^ C. □ 

We now introduce the following notion. Let 71 > be a constant. We shall 
refer to Rq as the following Robin function 

A x Ru(x,y) = -5(x -y) , infi\E, 
d v Rn(- 1 y)+~(iRn(;y) = , on 00, , (8.5) 

d„±R±(;y)- 1 ± (-)R±(;y)=0 , onS± , 

with j/ G f2 \ E. 

Lemma 8.2. Let R be the solution to (|4.1[) and let < r < ro, y £ SI \ E 6e 
suc/i i/iai B2r{y) C 0\S and dist(y,dfl) > ro- Then there exists a constant 
c r > depending on the a-priori data and on r only such that 

||-R(-,2/)IU°°(n\B r (i,)) < <V . (8.6) 

Proof. Let / 6 L~i~{VL) and let u 6 £T 1 (S1 \ E) be the weak solution to 

Au = f , in Q \ E , 

d^u + 71 u = , on dS7 , (8.7) 

d^u^ — 7 ± u ± = , on E . 

By Green's second formula the solution u can be represented as follows 

«(») = - / Rn{x,y)f(x)dx (8.8) 

where y 6 0\E. By the argument in Section 8.5] it follows that there exists 
a constant C > depending on the a-priori data only such that 

IMIl~ ( q\S) < C (jMI L2(n \S) + ll/ll L ^ M ) • ( 8 - 9 ) 

Moreover combining the weak formulation of problem (|8.7[) . the Poincare and 
the Holder inequalities we have that there exists a constant C > depending 
on the a-priori data only such that 

Furthermore, being 1 < ^ 2 we may infer that 

\\u\\ LH ^^C\\f\\ L ^ (8.11) 
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where C > is a constant depending on the a-priori data only. Hence inserting 
the above estimate in (|8.9|) we get that 



lL°°(n\s 



)<c\\s\\^ 



.12) 



where C > is a constant depending on the a-priori data only. 
Hence (U) and (f87T2|) yield to 



\\Rq(;V)\\ n±l _ 



sup 

11/11 a±l - = 1 



f(x)Rn(x,y)dx 



Finally by the weak Harnack inequality (see [22j Section 8.6]), we have that 
there exists a constant C r depending on the a-priori data and on r only such 
that 



\\Rn(;y)\\L^(n\B A y)) < C r \\R( Q ;y)\\ ^ 
Combining (|8.13[) and f|8. 13[) we obtain that 

\\Rn(-,y)\\L^(n\B r {y)) < Cr 



(Q\S) 



5.13) 



5.14) 



where C r > is a constant depending on the a-priori data only. Finally let us 
now consider the harmonic function uq(-) = R{-,y) — Rn(-,y) in 1 (S7 \ £). It 
follows that uq solves 

Au n = , in \ S , 

d u un = d v R(-,y) + JiRq(-, y) , on90 , (8.15) 
8 u ±Uq — 7 ± w^ = , on S ± . 

By (|8.14[) and by standard asymptotic estimate on the gradient of R(-,y) we 
get 



\\dvUn\\L°°(dn) < crj ™+7oC r 



(8.16) 



Classical estimates for harmonic functions leads to the existence of a constant 
C r > depending on the a-priori data only such that 



ll U n|| L oo( a \ s) < C r . 

Hence combining (|8.14[) and (|8.17[) we obtain the thesis. 



5.17) 

□ 



We now introduce the following notion. 

Let 70 be a constant. We shall denote with Rq the half space Robin function 



A x Rq(x, y) = -S(x - y) , in R» \ {x n = 0} 

d v ±Ro(-,y) -7 i?o = , on {x n = 0} , 



with y n ^0 and v+ = (0, ...,0,-1) and v~ = (0, ...,0, 1). 

Proposition 8.3. Let £ be a crack satisfying the a-priori assumption stated 
above. Let p > and let x 6 S p . Then there exists a constant C\,C2,C3 > 
depending on the a-priori data only such that 
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i) 



\V,R(z,y)\^co\z-y\ 1 - n , (8.19) 
for any y, z £ 1™. 

ii) 

\R( Zl y)-R Q ( Z ,y)\^^L\ Z -y\^+<* , (8.20) 



\V z R(z 7 y)-\7 z R (z,y)\^^\z-y\ 1 ~ n+a , (8.21) 

ro 

for any zeEfl B r (x) and for any y = hv(x) with < r < r~o, < h < f~o 
where fa — C3 min{ro, p} and 70 in (|8.18|) is such that 70 = 7 + (x). 

Proof. Without loss of generality we may assume that x = 0. Let po = 
\ min{ro, dist(x, and let $ e C 1 ' Q (B : aL,]R n ) be the map introduced in 

Theorem 13.21 In particular we have that for any < r < -^j it follows 

n+ n B d2 r(o) c *{B-{o)) cn+n B Bir (o) (8.22) 

where 9\ and 02 are the constants mentioned in Theorem 13.21 
We divide the proof in two steps. 

i) In the first step we shall prove that there exists a constant C\ > depending 
on the a-priori data only such that 

\V z R(z,y)\ ^dlz-j/l 1 -" for every z, y £ fi+ n B^g (0) , (8.23) 
the other cases being trivial. Let then z,y € f2+ fl Bg 2 m (0) and let £ = 

8M 

Denoting by 

A(C) = |det J D$(C)|( J D<I>- 1 )($(C))p$- 1 ) T ($(0) , (8.24) 
7+( C ) = 7 +($(0) , 7"(C) = 7"(*(0) (8-25) 
R(C, r,) = iZ(*(C), (8-26) 

it follows that 

Cil^l 2 < A(C)C • £ «S C 2 |^| 2 , V C € «), Ve e K" , (8.27) 

8M 



|A(Ci) - A(C 2 )| < C 3 |Ci - , VCi,C 2 eS^(0) (8.28) 

8M 

where Ci, Ci,Cz > are constants depending on the a-priori data only. Let us 
observe that R((, rf) satisfies 
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f -div c (A(C)V c J2(C,»7)) =S(C-v) , in 5^(0), 

< ~ ~ SM (8 29) 

\ A(C)V C i?(C,r?)V + 7+ (C)R(C,V) = 0, onB'™(0), V " ' 

where v' = (0,-^,0,1). 

Let v G H x (Q \ S) be the solution to the problem ([87Tj) and let «(C) = v(<E>(0). 
Since by Lemma T8. II we have that v(() ^ (7 in B7o_(0) then the quotient 

SM 

N(C,V) = ^^ (8-30) 

nC) 

is well defined there. 

Moreover, observing that A(£) = {cii j (C)) , i, J = 1, . . . , n is a symmetric matrix, 
we have that straightforward calculations lead to 



-div c (S(C)V c iV(C, v)) = S(C -V), in B™ (0) , 

_ 8M 

5(C)V c 7V(C,r?)-^' = 0, onB' (0) , 



5.31) 



where 5(0 = (MO) -J^ 2 (C)«,,i(C))- 
Writing for any C, £ _B_po_ (0) 



Ne(C,V) = 



N(C, v) if C = (0, Cn) is such that ( n ^ 
N((*,r)) if C = (0, CO is such that („ > 0, 

and for i = j with i, j = 1, . . . , n 



MO if C = (0, Cn) is such that Cn ^ 
bij(C) if C = (C'i Cn) is such that Cn > 0, 



whereas for i ^ j, i, j = 1 . . . , n we set 

MO = 



MO if C = (0, Cn) is such that Cn s$ 
-MC*) if C = (CCn) is such that Cn > 0, 

where C* = (C'i — Cn)- The first two are even and the third one is odd with 
respect to {Cn = 0}- In particular we have 

- div c (S e (C)V c iV e (C, r,)) = S(C -rj) + HC - V*) in (0) (8.32) 

where B%() = (6^.(0). 
Write 

N e (C,r)) = N e ((,v) + N e ((,V*) (8-33) 

where 

-div c ( J B e (C)V f iV e (C ) 7?)) = 5(C-^ in 5^(0). (8.34) 

Let r\ € B pq (0). By Lemma [8.21 we have that there exists a constant C > 
depending on the a-priori data only such that 

^••''Ww-n.*,^* - l8 - 35) 

V 128 128 
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We now consider the Green function G(C, r/) such that 

f -div f (S e (C)V C G(C,T?))=^(C-??) , inBrpoW 



G(C,r?)=0, ondil(O), (8 ' 36) 

8M 



with 77 € i? po (0). By the pointwise bound of G with the fundamental solution 
for the Laplace equation (see [29]) we infer that 

|G(C^)KG|C-^| 2 -" VCe-Br^O) , V /)eflpo(0) (8.37) 

8M IBM 

where G > is a constant depending on the a-priori data only. 
Let us define w((, rj) = N e ((, r/) — G((, 77), then we have 



dW c (B%O^MC,v)) = , in 5^(0) , 

8M 

w(C,r)) = N e ((,ri), on3B M (0) 



(8.38) 



Then by the bound in (|8.35[) and the maximum principle for solutions to equa- 
tions in divergence form we have that 

IH-^U^b m ^ C (8.39) 

8M 

where G > is a constant depending on the a-priori data only. 

Hence we may infer that there exists a constant G > depending on the a-priori 

data only such that 

\N e (Cv)\ <G|C-?y| 2 -" VCVqeflpn (0) , £^77. (8.40) 

Moreover recalling (|8.33[) we have that 

^(C^KGIC-*?! 2 -" VCVr/GB.go (0) , Z^v, (8-41) 

where C > is a constant depending on the a-priori data only. 

We observe that by Theorem 13. II we have that the function v(Q C in B~ P0 (0) 

where G > is a constant depending on the a-priori data only. Thus, by (|8.30[) 
and by (|8.41j) we get that 



\R{C,v)\<C\C-nr n VC,Vr ? eS^(0) , Z^r,, (8.42) 

where C > is a constant depending on the a-priori data only. 
Let h = dist(0, 77) = \r)\, then we have 



div c (A(C)V c i?(C,77)) = , in B£(0) 

A(C)V c i?(C,r?)-i/+7 + (C)%,r?)=0, onB'l(O) 



(8.43) 



By well-known regularity bounds for the Neumann problem (see for instance [21 
p. 667]) we have that 

l|V C i?(-,7?)|| LOO(B - (0)) < j\\R(;ri)\\L~ { B- m , (8-44) 
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where C > is a constant depending on the a-priori data only. By Theorem 
(regularity) we claim that there exists C G B~^(0) such that 

2 

\R(C,r,)\ = \\R(.,r,)\\ Lao — y (8.45) 
Then by (f8T42|) we find that 

^(■^^^jlC-Vl 2 -, (8-46) 

4 

where C > is a constant depending on the a-priori data only. 
On the other hand, noticing that 

\C-v\>\v\-\t\>h-± = ± (8.47) 

we obtain that 



IIW'^IU^) < CM 1 "" , (8.48) 

4 

where C > is a constant depending on the a-priori data only. 
Coming back to the original coordinates we have 



\V z R(z,y)\ = \D<f>-\z) T V c R($- 1 (z),<Z>- 1 (y))\^C 1 \z-y\ 1 - n , (8.49) 

where C\ > is a constant depending on the a-priori data only, 
ii) In the second step we shall achieve the desired asymptotic estimates. 
Arguing as in jTOl Proposition 3.4] we consider a function 8 £ C°°(M) such that 
< < 1, 0{t) = 1, for \t\ < 1, 6(t) = 0, for |i| > 2 and |^f| 2. Let us 
fix pi = min{j(8M )~ i , |} • and let us consider the following change of 
variables z = defined by 



C = 

Cn = *» - ^« (^) * (ff ) 



It can be verified that the map $ is a C 1,Q (Kl™, R") which satisfies the following 
properties 

§(Q-(0)) = Q+nQ pi (0) (8.50) 
c" 1 ^! - z 2 \ |* _1 (2i)-$ _1 (a;a)| < c|zi -z 2 | , V zi,z 2 6 M n , (8.51) 
li" 1 ^) - z\ «S — \z\ 1+a , V z £ R n , (8.52) 



\D$-\z)-I\< -^\z\ a , Vz£M" , (8.53) 
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where Q^O) = {( £ Q Pl (0) : ( n < 0} being Q Pl (0) the cube centered in O 
with sides of length 2pi and parallel to the coordinated axes and where c > 
is a constant depending on M and a only. 
Let us define the half cylinder C~ as 

C P1 ={zel": \z'\ < Pl , - Pl < z n < 0} . (8.54) 

For z, y G C~ , we have that R((, rf) = i?(<l>(£), $(ry)) is a solution to 



5.55) 



-div c (A(C)V c fl(C,77)2 = 5(C - 77) , in C" , 

A(C)V c i?(C, »j) V + 7+(C)i?(C, »?) = 0, on B' pi (0) , 

where z = "1(0,2/ = ^(??),7 + (C) = 7 + ($(0) an d where 

A(C) = |det J D$(C)|( J D$- 1 )(#(C))( J D*- 1 ) T ($(C))- 

Moreover, we observe that ^ is of class C a and A(0) = I. Let i?o(C: 7 7) be the 
fundamental solution introduced in (|8.18[) with 70 = 7 + (0). We notice that 
there exists a constant C > depending on the a-priori data only such that 
|7+(C) - 7o| < C\C\ a for any (' G fl^(0). 
Let us consider 

M(C,v) = R(Cv) - Ro(C,v) (8-56) 

which satisfies 

A C M(C, r?)) - div c ((J - A)(C)V CJ R(C, r?)) , in £7" , 

V c M(C,r/)- i / + 7o(C)M(C,r/)^_ 

= (I - A)Vi?(C, 7?) • 1/ + (70 - 7+(C))i?(C, »7) , on (0) . 

Let L > be such that C Bl(0). Thus by the representation formula over 
Bl(0) we get 

M(C,V) = I (/-l)(e)V £ %,C)V^ (e,r?)de + 

(7o - T+OW, 0), C)#o((£', 0), rj)d? + 



+ f (J-A)(ov c 5(^c)v e ieo(^i?)de + 

Ji3-(0)\C P - 

(^-i)(Ov f fi(e,o-^(e,»?) + 

a(B-(0))\B^(0) 

0„M(£, C)i?o(^, r/) - 9,i?o(e, V)M(C, 0M0 



>d(B-(0))\B' pi (0) 

For £|, \tj\ < ^- the last two integrals are bounded. Moreover, by (|8.23|) we have 
that 

\M((,v)\< C [ 1 + J c - \Z\ a \t-C\ 1 ~ n \t-*\ 1 ~ nd t) + 



c\ / \e\ a \(c,o)~c\ 2 - n \(e,o)-v\ 2 - n de 

V B ' n W , 
C(l + h + h + h + h) , 
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where C depends on the a-priori data only and 

h = I in-tf-^-vi 1 -^ 

•/c-n{|f|<4/>} 

h = [ me-cM^r^ 

^c-n{|f|>4fe} 



h= / \^\ a \(i'^)-cr n \{^Q)~vr n d£,' , 

JB> n (0)n{\e\<4h} 

h = ! ieri(^,o)-ci 2 " n i(e',o)- ?7 i 2 -"< , 

JB' P1 (o)n{\e\>ih} 

with h = |C — rj\. 

We bound I\ as follows 

h sC h a+2 - n [ |f~HC-£| 1 ""l'7-£| 1 ~ n d£< 

>/|fl<4 

< 4 a h a+2 - n [ iC-ll 1 ""^-!! 1 ""^ , (8-57) 

where £ = f , C = f ' ? ? = h • From standard bounds (see for instance [301 Chapter 
2]) it follows that 



/ lc-cT-"b7-eT""de<oo 

J|fl<4 



'lfl<4 

for any f , fj E M n , |C — ?} = 1. Thus we found that 

h «S c|C - ? /| Q+2 - n • 

Let us now consider li . We recall that by our hypothesis we have that |ry| = — r\ n . 
Let £ = (£', Cn) be such that < ^\Vn\- Then we have that h = |C — ^ \ \v\ 
from which we deduce that |£| ^ 2|£ — 77I and |£| ^ 4|£ — C\. Hence we obtain 
that 

h «S c / \£\ a+2 ~ 2n d£ < ch a+2 - n . 

J{\$\>4h} 

Treating analogously the integrals I3 and I4 we find that 

|M(C,t?)| ^C\C-v\ a+2 ~ n , (8-58) 

for any r\ = (0, • • • ,0, r) n ) such that < — rj n < ^ and for any |£| < ^ such 
that < \\rj n \ and where C > is a constant depending on the a-priori data 

only. 

Furthermore by Theorem l3.1l we have M(-, rf) £ C Q ( j_ ) where 1 = 

{C € C~ : I C| ^ |C«I ^ i 7 ?"}- Hence we can deduce that the above estimate 
remains valid for points \(\ ^ 4r such that £ = (C',0). 
We now go back to the original coordinates system. 
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Let z e $(B'n (0)) and let y = (0,y n ) with y„ g (-^,0), then since $ _1 (y) = y 

and since |$ _:L (2/)I = ~ $ _1 ( )l < - wc g ct b y (|83Tj) 

that 

c^\z\ < l*- 1 ^)! < ^(y) - ^(z)! + l^" 1 ^)! < c\y - z\. (8.59) 
On the other hand by (|8.52j) and by (|8.59[) we have that 

1$"^) - z\ < 4|z| 1+ « < - 2/| 1+Q • ( 8 - 6 °) 
Pa Pa 

Wc have that 

M(z,y) = - R (z,y) = M{$-\z), $" 1 (y)) + flo^" 1 ^), y) - . 

Then using (j8"3Tj) . (|53^1 . (j8~55)) and USSDH we find that 

|M(z,y)| < C|z-yr +2 -" + -^||Vi? (-,y)||^ (B , i(0)) |z-$- 1 (z)K 

< C|«- 1 ,| ot + 2 -» + ^-y| 1 -»|2 ! -y| 1 + a < 
Po 

*S ^-y| 2 +<*-«, (8.61) 

Po 

where C" > depends on the a-priori data only. 

We now estimate the gradient of M. Let z € ^(S'pi (0)) such that z = 3>(C) and 

4 

let /i = |C — y|. The following interpolation inequality holds 

||V c M0,y)|U~ (B , (0) ) <C\\M(.,y)\0 {BI ■ m \V ( M(-,v)\^ A0) , (8-62) 

4 ~4~ ~4~ 

where C > depends on the a-priori data only and 

|V c M(-,y)| a , fl / (0 )= sup . 

cc'^ (0) K-CI 

By the Holder continuity of V^-R and also of V^i?o we have that 
\VcM(;y)\ a ,B> h (o) < ( l|V c fl(-,2/)|| L oc (B / (0)) + ||V c i?o(-,y)||L-(s' i (o)) ) 

4 '*' \ 2 2 / 

sS Ch}- n - a (8.63) 

where C > depends on the a-priori data only. 
Hence combining (|8.61[) . (|8.62[) and (|8.63|) we get 

\V z M{<b-\z),y)\ < - y|( 2+Q -")(T^)| z _ y |(i-a-n)(T^) = 

Po +1 

= -§H*-l/| 1_n+3 & (8-64) 



Po +1 



33 



On the other hand we have 

\VM*~Hz),y)- V z R (z,y)\ < 

< \(D$-\z) T - I)VR (-,y)k- Hz) \ + \VR (-,y)k-i {z) - V z R (z,y)\ < 

< ^-\\VRo(;y)h- W o))\z - ^-\z)\ + \VRo(;v)\ a ,B> n (.o)\*-H*) - *l° < 
<-^\z- tf M* - + |* - y|- a+1 -"|* - y| (1+Q)Q < 

^-^\ z -y\^ n+a2 > ( 8 - 65 ) 

where C > depends on the a-priori data only. Thus by (|8.64[) and (|8.65|1 we 
obtain the thesis. □ 



Proposition 8.4. Let £ be a crack satisfying the a-priori assumption stated 
above. Let x £ f2 \ <9E and let y £ Q \ E such that \x — y\ ^ 2r$. Then there 
exist constants C, a > depending on the a-priori data only such that 

l V *^)l< g ( |^ r ~ 2+5 ) . ( 8 - 66 ) 
where r = j min{dist(x, <9£), \x — y\}. 

Proof. Let v £ H 1 (il\Y,) be the positive solution to (|8.ip introduced in Lemma 

R{x,y) 

18.11 We observe that N(x,y) = — '. solves 

chv x (u 2 (a:)V x jV(:r,?/)) = -<5(x - y) , infi\E , 
d„N(x,y)=0, onE. 

Let now a; be a point in fi \ 9E and let y £ f2 \ E. 

Let r be the radius defined as follows r = j min{dist(x, 9E), |x — y|}. Without 
loss of generality we may assume that x £ Let x' £ B r (x) n f2+ such that 
\x-x'\ §. 

By a change of variable argument we can deduce from Corollary 2.14 in [55] the 
following Holder continuity property of the Neumann function 

\N(x, y) - N(x', y)\ < C ^J^_ y ^ , (8.67) 

where C, S > depends on the a-priori data only. 
Moreover, being \x' — y\ ||x — y| we can deduce that 



|x - x' ld 

\x-y\ r 



\N(x,y) N(x',y)\ ^ C ,J" J 2+5 , (8.68) 



up to a possible replacing of the constant C in (|8.67[) . 

Let us now consider x £ dB r (x) (~1 fi+ and let x" £ B^{x) n fi+, then by the 
above estimate and by observing that \x — y\ ^ j\x — y\, we have that there 
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exists a constant C > depending on the a-priori data only such that 

\N{x",y)\ < \N(x",y)-N(x,y)\ + \N(x,y)\ ^ 

^ C ( \! X l y \S + ~ a +\*-y\ 2 - n ) . (8.69) 

Moreover, by the following local bound for the gradient, we have that there 
exists a constant C > depending on the a priori data only such that 

HVxJV(-,tf)|| L - ( fl s(lB) rin^ < ^\\N(-,y)\\ L ~ iB r( x) nTH) ( 8 - 7 °) 
By combining (|51)9"j) and (|57701) we find 

l|V^(,y)|U (Bi(x)n ^ < C ( - _^ 2+5 + ] -I± =g ) . (8.71) 

Next, being |x — y\ $J 2?*o we can find a constant C depending on the a-priori 
data only such that 



IIV.JVO.vJIL-^wrtQ < C lx _ yl n-2 + z ■ ( 8 - 72 ) 
Finally by the formal computation 

V x R(x, y) = \7 x N{x, y)v(x) + Vv{x)N{x, y) (8.73) 

and by analogous arguments of those applied above, the thesis follows. □ 

Let us consider O E Si HdV^ the point in (|4.14l) . We introduce a point O' G Si 
which is defined as follows by distinguishing two cases. 

• If O € Si is such that dist(0,9Ei) < then we consider a point 
O' € Si so that dist(0,C) = It follows that dist(C, <9£i) ^ % and 
dist(C,S 2 ) ^ 

• If O e Si is such that dist(0, <9£i) > %■ then we set O' = O. 
Proof of Proposition 14.61 We begin by recalling (|4.4p and we write 

|/(y,y)| > |S Sl (y,y)|-|S S2 (y,y)| . (8.74) 

First, we consider the term Ss 1 (y,y). We fix a radius p = min{ cc^, Coc^} and 
we observe that 
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(Si\s 2 )ns p (0') 



R+(-,y)d Ul Rt(;y)da 



\Ri {■,y)d Vl R 2 {-,y)da\ 



(Si\s 2 )ns p (0') 

|[iJi(-,tf)]i^i2a(-,»)|d<T 

(Ei\E 2 )\B„(0') 



(£i\£ 2 ) 



|[-R 2 (-,2/)9 l , 1 i?i(-,2/)]|dCT 



(8.75) 



Eins 2 



Let T{x, y) be the fundamental solution to the Laplace equation. Then we have 
that for any x £ B p (0') and for any y £ Br(0') 

f A x (R 2 (x,y)-T(x,y))=0 , in B^O') , 

where C is a constant depending on the a-priori data only. 
By the maximum principle for harmonic functions we get 

\R 2 (x,y)-T(x > y)\^Cp 2 - n V x £ B p (0'), V y £ Bz(0') . (8.76) 

By standard gradient estimates we have that 

\V x R 2 (x,y)-V x r(x,y)\^Cp 1 - n Vx£Be[0'), V y £ B § {0'). (8.77) 

Moreover observing that 



Ri(x, y) = Ro(x, y) + (Ri(x, y) - R (x, y)) 



we have that 



/si 



(s 1 \s 2 )nsp(0') 



Rf(x, y)9 v ( x )R 2 (x, y)da(x) 



(Ei\E 2 )nBp (C) 



Ro(x, y)d v ( x )T(x, y)da(x) 



\£ 2 )nB £ (0') 



\s 2 )nsp(0') 



\Ro(x,y)d u{x) (Rj(x,y) - T(x,y))\da(x) 
\(Rf(x,y) - R Q {x,y))d v ( x )T{x,y)\da{x) 



(8.78) 



V x R 2 (x, y) = V x T(x, y) + (V x i? 2 (x, y) - V x T(x, y)) (8.79) 



'(Ei\E 2 )nBp (0') 



\(Rf(x,y) - R {x,y))d u{x) (R%(x,y) - T(x,y))\da{x) 
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Let $ <E C 1 ' a (B_po_(0'),W l ) be the change of coordinates introduced in Theo- 

x 4M 

rem 13.21 then we have that = (0, . . . , 0, 1) + C?(|x'| Q ) and hence 



(Ei\2 2 )nB^(0') 



R (x, y)d u{x )T(x, y)da(x) 



-ca / - ?7| 2_2 ™|?T^' > Ci/i 2 -" - C 2 /i 2 -" +a 

'B' (*-i(0')) 

16M 

where ci, C2, Ci, C2 > are constants depending on the a-priori data only and 

We now consider the second term on the right hand side of (|8.75[) . We have 
that 

A x (R 1 (-,y))=0 , in B p (0') nO~, 

\Ri(;y)\\dB p (o')na- ^Cp 2 ~ n 

d u -R-(;y)- 7 -(-)R-(;y)=0 , in B p (0') n S~ . 
Hence by the weak maximum principle we have that 

2/) K<V" n in B p (0')nfl-, (8.80) 

where C > is a constant depending on the a priori data only. 
Then, by the asymptotic formulas (|8.19[) . (18.201) . we get 



/si 



(Si\s 2 )ns £ (0') 



'(£i\£ 2 )nB p (0') 
> C*i/i 2 "" - C* 2 /i 2 ~ n+Q 



#1 (%,y)dv( x )R 2 (x,y)\da(x) > 

C 3 |x-2/| 2 -"d^"da(x) 



(Si\s 2 )ns £ (0') 



(Si\s 2 )nSp(0') 
2 

C 5 



l i: \x-v\ J - n+a \x-y\ 1 - n d*(*) 



(Si\s 2 )nBp (O') «/i 



yja-n+adl-n^^ 



(8.81) 



(Si\s 2 )ns p (0') 



After straightforward calculation we observe that up to choosing the constant 
Co (in the definition of h and p) sufficiently small, we have that 



(Ei\E 2 )nB^(0') 



Rf (x, y)d v { x ) Rt (x, y)da(x] 



\Ri (x,y)d„( x )R 2 (x,y)\da(x) ^ ah 

'(Si\s 2 )nB p (0') 

where c\ > is a constant depending on the a-priori data only. 



2—n 
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We estimate the third term on the right hand side of (|8.75[) . By the asymptotic 
estimate (|8.19[) . we have that 

\R 1 ± (;y)\^C\x-y\ a - n ^C\^ 1 -h\' l ~ n on (Ei \ E 2 ) \ B p (0') 

where C > is a constant depending on the a priori data only. 
Moreover, by Proposition 18.41 we infer that, on (Si \ E2) \ B p (0') 

\& Vl R 2 {x, y)\ < C (dist(x, d^f^x - y\ 2 - n - 5 + \x- y^ . 

Hence by the integrability of dist(x,9E 2 ) over Si \ £ 2 we deduce that 

|[J2i(*,lO]ia„ (x )ife(a:,iO|tfcr(aO < C |t£ - hf~ 2n (8.82) 

(EAS 2 )\B p (0') 

where C > is a constant depending on the a priori data only. Finally, by 
(|8.19p and by the Robin boundary condition we get 



\R2{x,y)[d„ {x) R 1 (x,y)] 1 \da{x) + / \[R 2 (x,y)d Ul R 1 (x,y)]\dcr(x) 

2i\£ a JSinE 2 

^C\df -h\ 4 ~ 2n . 

Gathering together the above estimates we get 

\SvAy,v)\ > dh 2 - n - ci|^ - h\ 3 ~ 2n . (8.83) 

The upper bound 

\S Sa (v,y)\^&K -h\ 3 ~ 2n . (8.84) 

follows along the same lines of the arguments above. 

Combining the last two inequalities and (|8.74p we conclude the proof. 

□ 
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